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Abstract 

We calculate the total internal energy, total energy density and pressure, and the free energy 
for the neutrino and electromagnetic fields in Einstein and closed Priedmann cosmological models. 
The Casimir contributions to all these quantities are separated. The asymptotic expressions for 
both the total internal energy and free energy, and for the Casimir contributions to them are found 
in the limiting cases of low and high temperatures. It is shown that the neutrino field does not 
possess a classical limit at high temperature. As for the electromagnetic field, we demonstrate 
that the total internal energy has the classical contribution and the Casimir internal energy goes 
to the classical limit at high temperature. The respective Casimir free energy contains both linear 
and logarithmic terms with respect to the temperature. The total and Casimir entropies for the 
neutrino and electromagnetic fields at low temperature are also calculated and shown to be in 
agreement with the Nernst heat theorem. 

PACS numbers: 04.62. +v, 44.40. +a, 98.80.-k 
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I. INTRODUCTION 



The Casimir effect [l| manifests itself as some vacuum polarization energy and force 
due to material boundaries or nontrivial topology of space. This effect is presently well 
known owing to its multidisciplinary applications ranging from condensed matter physics 
and nanotechnology to quantum field theory, gravitation and cosmology (see monographs 
2|-|5|). During the last decade, numerous experiments on measuring the Casimir force have 
been performed as reported in reviews Q-a]- Their results were used to obtain stronger 
constraints on the parameters of the Yukawa-type corrections to Newtonian gravitational 
aw predicted in many extensions of the standard model (see the most recent results in Refs. 
9l-|l3|). Keeping in mind that measurements of the Casimir force are usually performed at 
room temperature, the problem of thermal Casimir effect became topical. In this regard it 
was found that the influence of free charge carriers in the boundary bodies on the thermal 
correction to the Casimir force is highly nontrivial and leads to complicated problems in 

li. 



quantum statistical physics js], ^ 

In cosmology, the Casimir effect arises not due to the material boundaries, but due to the 
identification conditions imposed on the wave functions by a nontrivial topology of space. 
Specifically, the Casimir energy density and pressure for the scalar field in the Einstein 
and closed Friedmann cosmological models with a topology x were found in Refs. 



15 



16| and in Ref. pJTjJ in the massless and massive case, respectively. The Casimir density 
and pressure for the neutrino and electromagnetic fields in the Einstein cosmological model 

n 

were derived in Ref. [18[. In succeeding years the Casimir effect was investigated in many 



topologically nontrivial spaces 



19N22|. This resulted in the development of new methods 



in mathematical physics and found prospective application in multidimensional cosmology 



23 



24| and in the problem of dark energy 



A considerable number of works was devoted to the thermal Casimir effect in cosmological 
models when not only zero-point photons, but thermal radiation as well was taken into 
account. Specifically, thermal Green's functions in Rindler, de Sitter and Schwarzschild 



spaces were considered in Ref. 



. In Refs. 



27 



28| the thermal stress-energy tensor of the 



scalar field in the Einstein cosmological model was investigated and the asymptotic behaviors 
at low and high temperature were found. Similar problem for a cosmological model with 
a 3-torus topology was solved in Ref. 29 1. In Ref. |3o| the results of Refs. 27, 28| were 
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reconsidered using another approach. It was shown that the total thermal stress-energy 
tensor of the scalar field in Einstein model contains the Casimir contribution which has the 
same asymptotic behavior at low and high temperature as in the case of two parallel ideal- 
metal plates. The total thermal stress-energy tensors of the neutrino and electromagnetic 
fields in Einstein cosmological model were considered in Ref . [Slj] . The expressions obtained 



were used 



31 



32| to determine the back reaction of the total thermal stress-energy tensor 



on the space-time by solving the Einstein equations, where this tensor would play the role 
of a source. 

It is common knowledge that at the very early stages the evolution of our universe has 
been going on at very high temperatures. Because of this, the studies of thermal Casimir 
effect in cosmological models mentioned above are of great physical significance. Keeping in 
mind, however, that massless fields of zero spin are not observed in nature, main attention in 
this respect should be paid to neutrino and electromagnetic fields. Considering that at large 
scales our universe is spatially homogeneous, here we derive exact expressions for the total 
thermal stress-energy tensor for these fields in Einstein and closed Friedmann cosmological 
models. We also find the total and Casimir free energy and the Casimir contributions to 
the stress-energy tensor (the latter were previously known only at zero temperature). We 
find the asymptotic behaviors of the obtained expressions in the limiting cases of low and 
high temperatures. For the neutrino case our results for the total energy and pressure are 
the same as in Ref. 3l|. For the electromagnetic case, however, there is a disagreement 
between our thermal stress-energy tensor and that obtained in Ref. SjJ- We explain this 
disagreement by the omission in jsil of a nontrivial contribution arising from the zeroth 
mode in the Poisson (or, equivalently, Abel-Plana) summation formula. The expression for 
the Casimir free energy and Casimir internal energy for the electromagnetic field obtained 
by us are in direct analogy with the familiar case of the thermal Casimir effect in an ideal 
metal spherical shell. For both the neutrino and electromagnetic fields we calculate the 
Casimir entropy and demonstrate that the third law of thermodynamics (the Nernst heat 
theorem) is followed. 

The structure of the paper is the following. In Sec. II we present general expressions for 
the free energy, internal energy and stress-energy tensor of the neutrino and electromagnetic 
fields in the Einstein cosmological model. Special attention is paid to the renormalization 
procedure and to the definitions of the Casimir contributions to the considered quantities 
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at nonzero temperature. Section III is devoted to the derivation of the total and Casimir 
free energy and internal energy for the neutrino field. The total and Casimir pressures are 
also considered. The asymptotic expressions of the obtained expressions at low and high 
temperatures are obtained and the validity of the Nernst heat theorem is verified. Similar 
results for the thermal Casimir effect in the electromagnetic case are presented in Sec. IV. 
Here, we concentrate on delicate points in the applicability of summation formulas and 
analogies with the case of electromagnetic field inside an ideal metal spherical shell. Our 
conclusions and discussions are presented in Sec. V. 

For convenience in comparisons with the classical limit (the case of high temperature 



33(1), below we retain the Planck constant h, the velocity of light c, and the Boltzmann 



constant ks in all equations. 



II. GENERAL EXPRESSIONS FOR THE TOTAL AND CASIMIR FREE ENER- 
GY AND INTERNAL ENERGY FOR NEUTRINO AND ELECTROMAGNETIC 
FIELDS 

We consider the four component neutrino field and the electromagnetic field in thermal 
equilibrium at some temperature T in the static Einstein cosmological model with a topology 
X R^. The metric of this model is given by 

ds^ = c^dr^ - a" [dr^ + sin^ r{de^ + sin^ Odip^)] , (1) 

where r, 9, Lp are dimensionless coordinates on a three-space of constant curvature +1 and r 
is the proper synchronous time. This is a spatially homogeneous model with a finite spatial 
volume V = 27r^a^, where a =const is the scale factor. All the results obtained below are 
applicable also in the conformally static closed Friedmann cosmological model with the scale 
factor depending on time (the concept of thermal equilibrium in nonstationary situations 



is discussed in Ref. 



3J]). In this case, however, there are additional contributions to the 



stress-energy tensor due to conformal anomaly and creation of particles [35|, |36|]. It has 
long been known that the three Friedmann cosmological models (open, quasi- Euclidean and 
closed) form the theoretical basis of modern astrophysics and cosmology. Because of this 
the results obtained below are not of only academic interest but can be used in cosmological 
applications. 
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The free energy of the neutrino (s = 1/2) and electromagnetic (s = 1) fields is given by 

f(^)(t) = eJ') + AF(^)(T), (2) 
where E^^^ is the zero-point energy at zero temperature and AF*^*)(T) is the thermal correc- 



tion. For the neutrino field we have 



18 



3l| 



(1/2) 



n=l 



1 + e 



[n + ^)c 



n=l 

For the electromagnetic field the respective expressions are the following 

oo 



18 



3l|: 



nc 



(3) 



(4) 



(5) 



n=l 



AF(^)(T) = 2kBTj2in' - l)ln [l - e-^^'^" V'^^^) 



n=l 



(6) 



Note that the zero-point energy E^^ given by Eqs. and (E]) is divergent, whereas the 
thermal correction AF^^^ (T) in Eqs. and is finite. It is notable also that the vacuum 
energy E^^'^^ is negative, as it should be in the spinor case. 

(s) 

The renormalization of the vacuum energy Eq is conventionally performed by subtracting 
the terms which are obtained from Eqs. (|3l) and ([5]) by replacing the discrete sums with 
integrations over a continuous variable 18|. As a result, the renormalized energies are 
defined by 



E, 



(1/2) 
O.ren 



-2h 



5^n(n+ 1)4^2) _ / 

n=l 



tit + l)uj\^'^^dt 



n=l 



(7) 



(8) 



The differences between the sums and t^ integrals can be simply evaluated by using the 
following Abel-Plana formulas 



\n] - I $ (t) dt 

n=l 



1 



(9) 



5 



n=l ^ 



$ (t) dt 



$ (it) - $ {-it) 



o2lTt 



+ 1 



(10) 



Equations and ( ITOi) are convenient in the cases of boson and fermion fields, respectively. 
The Abel-Plana formulas and fllOl) are valid for functions that are analytic in the right- 
half plane including the imaginary axis. 
The application of Eq. ffTOl) to ([7j) with 



$ ( n+ - 



$1 



(1/2) 



n -\ — 
2 



n{n + 1) I n + - 



1 

n -\ — 
2 



leads to 



18 



3l| 



(1/2) 
O.ren 



1 / 1 

4 I 2 



ir^c 



480a 



In a similar way, using Eq. (jH]) in ([8]) with 

$ (n) = <l>f ^ (n) = - 1) = 



n, 



one obtains |l8l. |31| 



E, 



(1) 

o.ren 



iihc 

120a 



(11) 



(12) 



(13) 



(14) 



As a result, the total free energy for the neutrino and electromagnetic fields in Einstein 
cosmological model is given by 



Fi:hT) = Ei^l^ + AE^^\T), 



(15) 



where all terms on the right-hand side are defined in Eqs. (jl]), ([6]), (fT2|) and ([T] 

Another important characteristic of the equilibrium state of the field at nonzero tempera- 
ture is the internal energy which is connected with the expectation value of the stress-energy 
tensor and the free energy in the following way 



ui:l{T) = v{Tt^') 

\ I tot 



d_ 
dT 



T 



(16) 



Using Eqs. @, (Q, f|T5|l and f lT6|) . we obtain for the internal energy of the neutrino and 
electromagnetic fields 



ui:l{T) = E^L+Au'~^\n 



(17) 



where 

,(1/2) 



A[/(V2)(T) =4;iV^^^^±ll^^, (18) 



n=l 



e 



A[/«(T) = 2hy ^^'"^^^"^ ■ (19) 



n=l 



Now we proceed with the definitions of the Casimir free energy and internal energy at 
nonzero temperature. These definitions go back to the Lifshitz theory |41] , where the Casimir 
free energy of the fiuctuating electromagnetic field between two semispaces separated by a 
gap of width a is obtained by subtracting the free energy of an unbounded Minkowski 
space. In other words, not only an (infinite) renormalization of the zero-point energy is 
done, but the thermal correction AF^^\T) undergoes finite renormalization by subtracting 
the contribution of the black-body radiation [s^. This contribution is proportional to the 
volume of the gap between semispaces. The same definition for the thermal Casimir effect 
was used in thermal quantum field theory for the configuration of two parallel ideal metal 
plates js, |4l|. What is more, for configurations with a finite volume, such as ideal metal 
rectangular boxes, it was shown 42| that to obtain the Casimir free energy one should 
subtract from AF^'\T) two more terms of quantum nature proportional to the surface 
area of the box and to the perimeter of its edges. This is in line with Ref. |[43|], which 
demonstrated that in the asymptotic limit of high temperature the total free energy in 
any restricted volume contains the following three types of terms depending on the Planck 
constant: 

is) {ksTY u) {ksTf u) {ksT? 

where the coefficients ot'Q \ oti \ 0-2'^ depend on the spin of the field and are expressed in 
terms of the heat kernel coefficients jsl. Some of these coefficients may be equal to zero. 



Specifically^ 



true that 



30| 



or a conformal massless scalar field in the Einstein cosmological model it is 



«r=-^^, «f)=o, 4°)=o. (21) 

It can be seen that in our case the terms (l20l) in the free energy are contained in the 
integral 

^^^\t)dt, (22) 



where the function ^'^\t) for the neutrino and electromagnetic fields in accordance with 
Eqs. (jl]) and ([6]) is given by 



<!>^i\t) = 2kBT{e - l)ln [1 - e-(^'"*/'^'=^^)] . (24) 



(23) 



Expanding the logarithms in power series and integrating, one obtains 44 1 

"90 {hcf ^ 12 he 



As a result, the Casimir free energy is defined as 

4^)(r) = 4J, + AFW(T), (27) 
where the Casimir thermal correction is given by 

poo 

AF^'\t) = AF^'\T) - ^^^'\t)dt. (28) 

Jo 

The definition fl28l) generalizes the approach previously followed in the Lifshitz theory and 
in thermal quantum field theory for ideal metal plates and rectangular boxes. As can be 
seen in Eq. ( l25l) and ( l26l) . for neutrino and electromagnetic fields in the Einstein model we 
have 

^(1/2) ^ J_!^y^ ^(1/2) ^ 0^ ^(1/2) _ ^ 



180 ' ' ' ^ 12 ' 

= «<" = 0, 4" = ^^. (29) 

In this manner, the definition of the thermal correction in the Casimir free energy presumes 
the subtraction of not only the contribution of the black-body radiation in free Minkowski 
space, as it is true for the scalar field [see Eq. but one more term of quantum nature 

which is present in the total free energy. In Sees. Ill and IV it is shown that the asymptotic 
expressions of the total free energy F^^^ at high T do not contain any more power-type terms 
of quantum nature. 
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The Casimir contribution to the internal energy (and respective Casimir contribution to 
the thermal stress-energy tensor) can be defined in a similar way. It is given by 

4^)(T) = i?g,„ + Af/(f)(T), (30) 

where 

/•oo 

AU^c\t) = AU^'\T) - ^i'\t)dt. (31) 

Jo 

In accordance with Eqs. (fTSll and (fT9|) . the function ^l^\t) for the neutrino and electromag- 
netic fields is given by 

, 2hc t(t^ - 1) 
^3^^(^) = — e(^-cW)-l - (33) 



The integrals subtracted in Eq. fl3T]) are simply calculated as 

TT^a^ (A;bT)4 TT^a (/ 
^0 {hcf 12 ~ he 





It is easily seen that the integrals ( l25l) and ( |26|) subtracted from the total free energy and 
the respective integrals ( l33l) and ( IMl) subtracted from the total internal energy satisfy the 
same Eq. (fT6l) as the total quantities. 

In the end of this section, we consider other components of the stress-energy tensor diffe- 
rent from 00-component defined in Eq. (1161) . In the spatially homogeneous isotropic metrics 
of the space-time under consideration the stress-energy tensor is diagonal and has equal 
spatial components 

pi:hT) = - (rt^') , (36) 



tot 

where (T) is the pressure. Both the total and Casimir pressures are expressed using the 
respective free energies as 

= p'cHt) = (37) 

9 



Considering Eqs. (H]), (E]), (fT2|) . ( !T^ . and (fT5|) . we obtain for the total pressure defined 
in Eqs. fl37|) . the following results 



(1/2) . . ^ 2/1 >^ n(n + 



Act I J J 7207r2a4 ^ 37r2a3 ^ .(a^'/fcaT) _ i ' ^''''^ 

n=l 

Comparing Eq. (EH]) with Eqs. ([I6])-([I8]) and Eq. ^ with Eqs. ([16]), ([I7D and ([19]), we 
get the equation of state 

= \e[i{n (40) 
where the total energy density is just given by 

.S(r) = (rr) (41) 

\ /tot 1/ 

in accordance with Eq. f[T6l) . Using Eqs. f[25l) . f[26]) . f[34l) and f[35]) . it can be easily verified 
that the same equation of state is satisfied for the Casimir quantities defined in Eqs. ([271) . 
(I2H]), dsn]), dSI]) and dSB: 

p(f)(T) = \e^{T). (42) 

Here, the Casimir energy density is the 00-component of the diagonal Casimir stress-energy 
tensor 



4'l(r) = (t<"»)^ ^ (43) 



The other components of this tensor are given by —P^ iT). 

III. CALCULATION OF THE TOTAL AND CASIMIR FREE ENERGY AND IN- 
TERNAL ENERGY FOR NEUTRINO FIELD 

Here, we calculate the free energy and internal energy for the neutrino field in Einstein 
and Friedmann cosmological models. We begin from the representation of the Casimir free 
energy Fq^'^\t) as given by Eq. ([271) with s = 1/2, where the energy at zero temperature 
is expressed by Eq. ([T2|) and the Casimir thermal correction by Eqs. ([1]), ([23l) and ([28l) . 
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It is convenient to calculate the difference between the sum in Eq. (jl]) and the integral in 
Eq. ([28]) by using the Abel-Plana formula ([lO]) with <l>(t) = ^2^^\t) defined in Eq. ([23]). 
This function is analytic in the right-half plane including the imaginary frequency axis. 
Expanding the logarithm in power series, one obtains 



(t) = -AksT --)f2 ( ^)"'"' e-(r»cniMBT) (44) 



n=l 

and 



n=l m=l ^ 



_ ^ _8,,^j, A. tip. (45) 

^ ' ?1=1 

Then, calculating the integral on the right-hand side of Eq. (fTOj) . we arrive at the represen- 
tation of the Casimir free energy in terms of the double sum: 

{hcn/akBTf - 3(27rm)2 
[{hcn/akBTy + {2nmy]^ 

8 {hcn/aksTy + (27rm)2 j ' ^^^^ 
In Eq. ( H6|) one can carry out the summation either in m or in n. The resulting repre- 
sentations for the Casimir free energy are, strictly speaking, equivalent, but convenient for 
obtaining the asymptotic limits in the case of low and high temperature, respectively. The 
summation in m can be performed by using the formula 

x'^ — Sy'^m'^ 1 



1 



m=l 



^2 _|_ ^2^2)3 g(^jr2 _|_ y2^2) 

y — nxcschlrcx/y) 



{4 - 7r^a;^[3 + cosh(27rx/2/)]csch^(7rx/y)} 



(47) 



IGx'^y 

where x = hcn/akBT and y = 27T. By performing the summation in m in Eq. ( I46l) with the 

help of Eq. fH7|) , one arrives at 

^(1/2) I7hc ^ 7n^a' {ksT)^ n^a {ksTf 
= 7^ + 



480a 90 {hcf 12 he 



^ ^ n sinh^(ncn/2aA;BT)' ^ ' 

From this equation, in the low temperature limit, where hc/akBT — ^ oo, we have 

^ ^ ' 480a 90 {hey 12 /ic 

- 8A;BTe-(3'^^/2afcsT)_ (4g) 
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This result is similar to the low-temperature behavior of the Casimir free energy in the case 
of two parallel ideal metal plates 

From Eq. (148 p it is easy to obtain the exact expression for the total free energy. Using 
Eqs. f lT5]) . fl2^ and fl2S]) one has all power- type contributions in T canceled, and then the 
total free energy is given by 

17hc , 1 



f(1/2) 



(T) 



480a 



^ n 



(50) 



sinh {hcn/2akBT) 

Thus, the thermal correction in the total free energy at low temperature is exponentially 
small: 



.(1/2) 



tot 



(T) 



480a 



8A;Rre-(=^'^^/2afeBT)_ 



(51) 



To obtain the asymptotic behavior of the Casimir and total free energy at high tem- 
perature, we perform the summation in n in Eq. ( I46l) first. This leads to the following 
representation for the Casimir free energy: 



7(1/2) 
c 



(T): 

{hcf 



17 he {kBTf 



+ 



E 



\m+l 



480a 167r4 {hcf ^_^m^smh\27c^makBT/hc 

Air'^makBT 



+3 



(akBT)^ 
+47r^m 
{he) 



,42, 

(6 + TT^m^) sinh 



(2 + 7r'^m^) 



2'K'^makBT 



cosh 



he 



he 



^^^^ (2 + TT^m^) + 47r^ .^'^„ m sinh 



(aA;Br)^ 



akBT 



An'^makBT 
he 



(6 + vr m ) sinh 



Gn^makBT 



{akBTY ' ' \ he 
The same expression can be equivalently rewritten in the form 



(52) 



7(1/2) 
c 



17he {kBTf 



\m+l 



m=l 



X 



X 



480a ^47r4 {he)^ 

{heY , 9 9x 9 , / 2'iT^makBT 

^ ' r(6 + 7r2m2) + 2m7r2csch ' 



{akBTf 
, hem ( hem 



he 



IT 



+ 4coth 



4m^7r'^csch 



akBT \akBT 

2Ti'^makBT 



2Tx'^makBT 
he 



+ 2 



{he) 



{akBTf 



+ 2mV^ 



he 

From Eqs. (l52l) or (!53i) in the limit of high temperature one obtains 



2 2(^5^)' -( 2n akgT/hc) 



{hef 



(53) 



(54) 
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Note that in the hmit of high T the zero-temperature contribution to the Casimir free energy 
was canceled by the corresponding term in the thermal correction. We emphasize that the 
Casimir free energy at high T is exponentially small. It does not contain the classical term 
proportional to fc^T (see Section IV), as is expected for a spinor field. 

The total free energy is obtained by adding Eq. to Eqs. (13^ or (1351) . The expres- 
sions obtained for f/o/^^(T) are equivalent to the more compact Eq. (!50|) . The asymptotic 
behavior of the total free energy at high T is obtained by adding Eqs. fl25|) and fl5^ . 

A quantity closely related to the free energy is the entropy. The total and the Casimir 
entropies for both the neutrino and electromagnetic fields are defined as 



(55) 



The exact expression for the Casimir entropy of the neutrino field can be obtained from Eq. 
(l48l) . It is given by 



S^c'^\t) 



2, akBT 

TT ks 

She 



1 _ Mvr^ {akBT f 

2 ^ ^5 {hey 



n=l 
,4 



n 



csclr^ {hen /2akBT) + 



3hen 
2akBT 



cosh 



hen 
2akBT 



X csch {hen/2akBT 



(56) 



When the temperature vanishes, the asymptotic expression for the Casimir entropy is given 
by 



+ 



TC^kB- 

12he 
aT 



She 



1 _ Utt^ {akBTf 

2 ~ 



15 {he) 



-{3hc/2akBT) 



(57) 



As is seen from Eq. flFTj) . the Casimir entropy goes to zero when T goes to zero, i.e., the third 



law of thermodynamics (the Nernst heat theorem) is satisfied 39|, 4^. The total entropy 
is obtained from Eqs. (!50|) and (!55l) . It is given by Eq. (!56l) with the first term on the 
right-hand side omitted. The asymptotic behavior of the total entropy at low T is given by 
the exponentially small term in Eq. ( 157|) . Thus, the Nernst heat theorem also holds for the 
total entropy of neutrino field. 

Now we turn to the consideration of the internal energy. The total internal energy for the 
neutrino field is given by Eqs. ( IT71) and (fT8|) . The Casimir internal energy can be calculated 



13 



by Eqs. fl30|l and fl3T]) with the function ^^^'^\t) defined in Eq. fl32|) . The calculation follows 
the same steps as shown above for the free energy and uses the Abel-Plana formula (fTO|) . 
The same results can also be obtained from the respective expressions for the free energy by 
using Eq. (fT6l) and similar equations for the Casimir quantities. Thus, from Eqs. (1501) and 
( |T6|) . we find the following result for the total internal energy: 

^ - 480a + 2a smh\hcn/2akBTy ^ ^ 

In the limit of low temperature, from Eq. fl55]) . one obtains 



In the limit of high temperature the asymptotic behavior of the total internal energy is 



30 {hcf 12 he 



4 



Dividing both sides of Eqs. (l59l) and (l60l) by the spatial volume V , we get 



-t'^(^) = ^i7^ + ^-'^"^'"^^^^ (61) 
9o07r^a^ vr^a^ 



Tn^kBTY 1 (A;bT)^ 



60 (nc)3 24a2 he 

^^^2{±^^-i2.^aksT/nc) (g2) 

at low and high temperature, respectively. These coincide with the asymptotic behavior 



of the energy density for the neutrino field found in Ref. 3ll] (note that the exponentially 



small terms were omitted in 31|). 

The exact expression for the Casimir internal energy is obtained by subtracting Eq. ( 134|) 
from Eq. fl58l) . As a result, at low and high temperature we have, respectively, 

(1/2). . ^ 11^ _ 77r^a' [ksTY -n'^a {ksTf 
^ ^ ^ 480a 30 {hef 12 he 

^l2hc^-isnc/2ak,T)^ (63) 
a 

U^^"\t) = 87r^a3%^e-(2-'"'=-^/^''=). (64) 
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Thus, for the neutrino field the Casimir internal energy at high T is exponentially small. 
The limiting values for the Casimir energy density are obtained from Eqs. fl63l) and f lM|) 
after the division by V . As to the total and Casimir pressures, they can be obtained from 
the equations of state fllUl) and (jUj). 



IV. CALCULATION OF THE FREE ENERGY AND INTERNAL ENERGY FOR 
ELECTROMAGNETIC FIELD 

The case of the electromagnetic field is more complicated because the function $2^'* (t) defined 
in Eq. ( 12^ goes to infinity when t goes to zero. This prevents direct application of the Abel- 
Plana formula ([9]) for the calculation of the Casimir free energy of the electromagnetic field 
defined in Eqs. ( 127|) and (l28l) . A similar problem arises in the calculation of the total and 
Casimir internal energies given by Eqs. ffT7|) . f lT^ and f l5U]) . fl^ . respectively. The point is 
that the function (^^^\t) in Eq. f l33p . determining the thermal correction f ll9p in the internal 
energy, has poles along the imaginary frequency axis. This also prevents the application of 
the Abel-Plana formula in its simplest form ([9]). 

We start from the total free energy ([2]), ([6]), and expand the logarithm into the power 
series to obtain 

^^^^^^ = - 2^^^ S ^ - l)e-^^^""^/"'=-^^ (65) 

m=l 71=1 

Here, the sum with respect to n can be calculated leading to 

*°* ^ ' 120a ^ ^(^hcm/aksT _ ^ 3 ^ ' 

m=l ^ ' 

In accordance with Eq. fl251) the Casimir free energy for the electromagnetic field is obtained 
subtraction of Eq. fl26l) from Eq. fl66|) : 

^ ^ ' 120a 45 {hcf 3 /ic 

-2fcBTy , . , , ^ — — r. (67) 

Equations (|66!) and fl671) are convenient for obtaining the asymptotic expressions at low 
temperature. For the total free energy one obtains from Eq. (!66|l 



FtS (^) = ^ - ^ksTe-^-^^l^'-^^ . (68) 
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It is seen that total correction to the energy at zero temperature is exponentially small. 
The low temperature behavior of the Casimir free energy is given by the difference between 
Eqs. (IBS]) and In this case, there are power- type corrections to the energy at zero 

temperature. 

Before considering the asymptotic behavior of the free energy at high T, we obtain the 
exact expression for the internal energy and its asymptotic behaviors at low and high tem- 
peratures. Using the expansion in power series in Eq. f fTOj) . we can rewrite Eq. f|T7|) in the 
form 



oo oo 



uilhT) = ^ + — E - l)e-^'--/^'-^\ (69) 



120a a 

m=l n=l 



By performing the summation in n we obtain the following exact expression: 

120a a (e^cm/aksT _ i^A- y'") 

m=l ^ ' 



The Casimir internal energy is given by Eqs. (130|) and (13T]) . Taking into account Eq. ( 135|) 
we have 

(1) ^ 11^ _ 27rV {k^Tf n^iksTf 
^' ^ ' 120a 15 {hey 3 he 

I2hc ^2hcm/akBT 

^ ^ Uhcm/akBT _ ]^-\4" ^'^^^ 

m=l ^ ' 

Equations ( 170]) and (17T]) are convenient for obtaining the asymptotic expressions at low T . 
Thus, from Eq. (!70!) it follows 

After subtraction of Eq. (!35l) from Eq. (1721) . the asymptotic expression for the Casimir 
internal energy at low T is obtained. All these expressions for the internal energy are 
connected with respective expressions for the free energy, obtained above, by Eq. (fT6|) and 
a similar equation for the Casimir quantities. 

Now we turn to the Casimir internal energy in the form of Eqs. (1301) and (I3T1) . and 
calculate it in a different way, which allows us to obtain the high temperature limit. The poles 
of the function $3^''(t) defined in Eq. (|33|) are located at the points t = iti = 27rilakBT/hc, 
where / = ±1,±2,... . In addition, at t = this function takes a nonzero value $^^^(0) = 
—2kBT. The presence of the poles along the imaginary frequency axis makes the standard 
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Abel-Plana formula inapplicable. There is, however, a generalization of this formula for the 



case when there are poles on the imaginary frequency axis. It reads 



37| 



Hn) - / mdt = --m Res 

n=l ^ I ' 



sin(7rt) 



+1 I i . — -dt. 

'o 



- 1 



(73) 



Here, it is assumed that summation is done only over the positive / and the function 
satisfies the condition 

1 



= + o 



when t iti. 



(74) 



We apply Eq. ( I73|) to the function $(t) = ^^p{t) defined in Eq. ( 133|) . In doing so, we 



take into account that 



TT y Res 

^ t=it. 



1=1 



sin(7rt) 



(he) 



1=1 



1=1 



he 



and 



2 he 



it{t^ + 1). 



(75) 



(76) 



Calculating the integral on the right-hand side of Eq. fl73|) . 



^^^\it)-^t\-itl^^^_ ll he 



(77) 

.(1) 



,0 e^^^* - 1 120 a 

we find that it cancels the contributions of the Casimir energy at zero temperature Eq^^^ 
defined in Eq. f lT4|) . Then from Eqs. f l30|) and fl3T|) one obtains 



(nc)3 



^Aw^lakBT/hc _ 



I 



he 



(78) 



This is an alternative exact expression for the Casimir internal energy. Numerical compu- 
tations show that is gives the same values of UqI^{T) as Eq. (ITTj) over the entire range of 
temperatures. The total internal energy u[\^{T) is obtained by adding Eq. (|35|) to Eq. (!78|) . 
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Equation (178|) is convenient for obtaining the asymptotic limit at high temperature. Thus, 
from this equation the asymptotic expression for the Casimir internal energy at high T is 
given by 



u'^^\T) = kBT + 



32n^a%kBTy 



(79) 



The asymptotic expression at high T for the total internal energy is obtained from Eq. ([79 
by adding Eq. ( 135|1 : 



15 {he) 



he 



+ kBT 



I 327r a^kBT) (4^2 ^^^y/M 
^ (^c)3 



(80) 



From Eq. f l79p it is seen that the electromagnetic Casimir internal energy at high tempera- 



ture has a classical limit 



33| . which does not depend on h and c. The total internal energy 



f lHOj) also contains a classical term equal to kBT. The respective asymptotic expression for 
the total internal energy density is 

uilkT) -K^kBTY 1 {kBTf 



-(1) 



(T) 



kBT 



V 



15 {he) 



Qa? he 



27r2a3 



{hef 



(81) 



Note that the high temperature behavior of elol obtained in Ref. |3l| does not contain the 
classical term. This is explained by the omission of the term $3^^(0)/2 = -fc^T in the 
Poisson summation formula (Eq. (4) of Ref. 3l|), which is equivalent to the Abel-Plana 
formula used here. 

The obtained expression f l79p for the Casimir internal energy at high T can be used for 
the determination of the high temperature behavior of the Casimir free energy, which is still 
unknown. Using the thermodynamic connection (|T6|) between the free energy and internal 
energy, we have 

U^c\t) 



F^c\t) = -T 



J^2 



-dT + kBR 



(82) 



where R is an arbitrary dimensionless constant independent of T. Substituting Eq. (17^ in 
Eq. (182!) and keeping only the high order terms, one arrives at 



F^^\t) = -kBTln^^^ - RkBT 



akBT 



+87r"a 



he 

2 2ikBT)' 



{hey 



(83) 
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The value of the constant R can be determined from the exact expression for F^^\t) see 
Eq. (!67|) ]. Computations show that the exact expression (!67|) leads to the same values of 
the Casimir free energy as the asymptotic expression fl83|) up to six significant figures for 
the values of parameter aksT/hc > 1 when R = 1.77698. Thus, the Casimir free energy 
of the electromagnetic field at low T contains not only the classical (entropic) term, but 
the logarithmic contribution as well. The same holds for the Casimir free energy inside an 
ideal metal spherical shell js], Q, 46 1 (at zero temperature the Casimir ener gy for fields of 



different spins in spherically symmetric cavities was considered in Refs. |47|, |48|). The low 
temperature behavior of the total free energy is obtained from Eq. ( |83l) by adding Eq. (|26l) . 

In the end of this section we consider the total and Casimir entropy for the electromagnetic 
field in Einstein and closed Friedmann cosmological model. At low temperature from Eqs. 
(!55|) and (!68|) one obtains the following main contribution to the total entropy: 

SilliT) = l^e-(^'^^M-^). (84) 

For the Casimir entropy at low T using Eq. ( 126|) we arrive at 



4'^(T) =2'K''kB 



aksT 



She 

12hc 



15 {hcf 



^-_^^-i2Kclak^T) _ (85) 

aT 

As is seen from Eqs. f lH^ and fl85l) . both the total and Casimir entropies go to zero when 
the temperature vanishes, in accordance with the Nernst heat theorem. 

Using the above results, the expressions for the total and Casimir pressures of electro- 
magnetic field and their behaviors at low and high temperature can be obtained from Eqs. 
gOD and gID. 



V. CONCLUSIONS AND DISCUSSION 



In the foregoing we have investigated the thermal Casimir effect for the neutrino and 
electromagnetic fields in the Einstein cosmological model. The results obtained are also 
valid in the closed Friedmann cosmological model where they should be considered as com- 
plementary to the terms describing the vacuum polarization and particle creation caused 
by the nonstationary regime of the metric. For both fields under consideration we found 
general expressions for the total internal energy, energy density and pressure (considered 
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earlier) and for the total free energy. In all cases, we separated the Casimir contributions 
from the obtained quantities by means of an additional subtraction procedure. 

Both the total and Casimir internal energy and free energy were represented in the form 
of single sums. The asymptotic expressions for these sums were found in the limiting cases 
of low and high temperature. For the neutrino field, our results for the total internal energy 
(energy density) are in agreement with those obtained in Ref. jsi]. We have also found 



3l| . Our results for the Casimir 



the exponentially small corrections to the results of Ref. 
free energy of neutrino field are similar to those obtained for configurations with material 
boundaries. Specifically, we have shown that the Casimir free energy of neutrino field does 
not possess a classical limit at high temperature and the Casimir entropy satisfies the Nernst 
heat theorem. 

For the electromagnetic field, our calculation result for the total internal energy differs 
from the result of Ref. 3l| by a classical term which arises from the contribution of zero 



arguments in the Poisson and Abel-Plana summation formulas. We have shown that for 
the electromagnetic field at high temperature the total internal energy includes not only the 
Planck-type terms, as it was believed previously, but also the linear in temperature classical 
term. For the Casimir free energy of the electromagnetic field in Einstein cosmological model 
we have proved that there are both classical and logarithmic in temperature terms in the 
limit of high temperature, as it holds for the thermal Casimir effect inside an ideal metal 
sphere. Both the total and Casimir entropies were demonstrated to be in agreement with 
the Nernst heat theorem. 

In future, it would be interesting to determine the influence of the revealed classical term 
on the cosmological evolution, where the total internal energy plays the role of a source, and 
consider mult i- dimensional generalizations of the obtained results for application in brane- 



world scenarios. Specifically, it seems advantageous to generalize the results of Refs. 23l. |24| 
to the case of nonzero temperature. 
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